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§5.2 CONSTITUTIVE EQUATIONS 


§5.1. Introduction 


Recall from the previous Lecture the following connections between various quantities that appear 
in continuum structural mechanics: 


; MP ; 

internal forces > stresses => strains = displacements => size & shape changes (5.1) 
... MP 
displacements = strains = stresses = internal forces (5.2) 


Of these, we have studied mechanical stresses in Lecture 1 and strains in Lecture 4. How are they 
linked? Through the material properties of the structural body. This is pictured by the ‘MP’ symbol 
above the appropriate arrow connectors. Material behavior is mathematically characterized by the 
so-called constitutive equations, also called material laws. 


§5.2. Constitutive Equations 


In this Lecture we will restrict detailed examination of constitutive behavior to elastic isotropic 
materials. More complex behavior (for example: orthotropy, plasticity, viscoelasticity, and fracture) 
are studied in senior and graduate level courses in Structural and Solid mechanics. 


§5.2.1. Material Behavior Assumptions 


There is a very wide range of materials used for structures, with drastically different behavior. In 
addition the same material can go through different response regimes: elastic, plastic, viscoelastic, 
cracking and localization, fracture. As noted above, we will restrict our attention to a very specific 
material class and response regime by making the following behavioral assumptions. 


1. Macroscopic Model. The material is mathematically modeled as a continuum body. Features 
at the meso, micro and nano levels: crystal grains, molecules, and atoms, are ignored. 


2. Elasticity. This means the stress-strain response is reversible and consequently the material 
has a preferred natural state. This state is assumed to be taken in the absence of loads at a 
reference temperature. By convention we will say that the material is then unstressed and 
undeformed. On applying loads, and possibly temperature changes, the material develops 
nonzero stresses and strains, and moves to occupy a deformed configuration. 


3. Linearity. The relationship between strains and stresses is linear. Doubling stresses doubles 
strains, and viceversa. 


4. Isotropy. The properties of the material are independent of direction. This is a good assumption 
for materials such as metals, concrete, plastics, etc. Itis not adequate for heterogenous mixtures 
such as composites or reinforced concrete, which are anisotropic by nature. The substantial 
complications introduced by anisotropic behavior justifies its exclusion from an introductory 
treatment. 


5. Small Strains. Deformations are considered so small that changes of geometry are neglected 
as the loads are applied. Violation of this assumption requires the introduction of nonlinear 
relations between displacements and strains. This is necessary for highly deformable materials 
such as rubber (more generally, polymers). Inclusion of nonlinear behavior significantly 
complicates the constitutive equations and is therefore left for advanced courses. 
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FIGURE 5.1. Typical tension test behavior of mild steel, which displays a well defined yield 
point and extensive yield region. 
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FIGURE 5.2. Three material response “flavors” as displayed in a tension test. 
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FIGURE 5.3. Different steel grades have approximately the same elastic modulus, but very 
different post-elastic behavior. 
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§5.3 CHARACTERIZING A LINEARLY ELASTIC ISOTROPIC MATERIAL 


$5.2.2. The Tension Test Revisited 


The first acquaintance of an engineering student with lab-controlled material behavior is usually 
through tension tests carried out during the first Mechanics, Statics and Structures sophomore 
course. Test results are usually displayed as axial nominal strain versus axial nominal strain, as 
illustrated in Figure 5.1 for a mild steel specimen taken up to failure. Several response regions 
are indicated there: linearly elastic, yield, strain hardening, localization and failure. These are 
discussed in the aforementioned course, and studied further in courses on Aerospace Materials. 
It is sufficient to note here that we shall be mostly concerned with the linearly elastic region that 
occurs before yield. In that region the one-dimensional (1D) Hooke’s law is assumed to hold. 


Material behavior may depart significantly from that shown in Figure 5.1. Three distinct flavors: 
brittle, moderately ductile and nonlinear-from-start, are shown schematically in Figure 5.2 . Brittle 
materials such as glass, rock, ceramics, concrete-under-tension, etc., exhibit primarily linear be- 
havior up to near failure by fracture. Metallic alloys used in aerospace, such as Aluminum and 
Titanium alloys, display moderately ductile behavior, without a well defined yield point and yield 
region: the stress-strain curve gradually turns down finally dropping to failure. Some materials, 
such as rubbers and polymers, exhibit strong nonlinear behavior from the start. Although such 
materials may be elastic there is no easily identifiable linearly elastic region. 


Even for a well known material such as steel, the tension test behavior can vary significantly 
depending on combination with other components. Figure 5.3 sketches the response of mild steel 
with high-strength steel used in critical structural components, and with tool steel. Mild steel is 
highly ductile and clearly exhibits an extensive yield region. Hi-strength steel is less ductile and 
does not show a well defined yield point. Tool steel has little ductility, and its behavior displays 
features associated with brittle materials. The trade off between ductility and strength is typical. 
Note, however, that all three grades of steel have approximately the same elastic modulus, which is 
the slope of the stress-strain line in the linear region of the tension test. 


§5.3. Characterizing a Linearly Elastic Isotropic Material 


For an isotropic material in the linearly elastic region of its response, four numerical properties are 
sufficient to establish constitutive equations. Those equations are associated with the well known 
Hooke’s law, originally enunciated by Robert Hooke by 1660 for a spring, and later expressed 
in terms of stresses and strains once those concepts appeared in the XIX Century. These four 
properties: E, v, G and q, are tabulated in Figure 5.4. 


$5.3.1. Determination Of Elastic Modulus and Poisson’s Ratio 


The experimental determination of the elastic modulus E and Poisson’s ratio v makes use of a 
uniaxial tension test specimen such as the one pictured in Figure 5.5. See slides for operational 
details. 


§5.3.2. Determination Of Shear Modulus 


The experimental determination of the shear modulus G makes use of a torsion test specimen such 
as the one pictured in Figure 5.6. See slides for operational details. 
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E _ Elastic modulus, a.k.a. Young's modulus 
Physical dimension: stress=force/area (e.g. ksi) 


Vv Poisson's ratio 
Physical dimension: dimensionless (just a number) 


G_ Shear modulus, a.k.a. modulus of rigidity 
Physical dimension: stress=force/area (e.g. MPa) 


a Coefficient of thermal expansion 
Physical dimension: 1/degree (e.g., 1/°C) 


E, v and G are not independent. They are linked by 


E=2G(14v), G=E/Q204\)), v=E/2G)-1 


FIGURE 5.4. Four properties that fully characterize the thermomechanical response 
of an isotropic material in the linearly elastic range. 


§5.3.3. Determination Of Thermal Expansion Coefficient 


The experimental determination of the thermal expansion coefficient a can be made byt using a 
uniaxial tension test specimen such as the one pictured in Figure 5.7. See slides for operational 
details. 


§5.4. Hooke’s Law in 1D 


Once the values of E, v, G anda are experimentally determined (for widely used structural materials 
they can be simply read off a manual), they can be used to construct thermoelastic constitutive 
equations that link stresses and strains as described in the following subsections. 


§5.4.1. Elastic Modulus And Poisson’s Ratio In 1D Stress State 


The one-dimensional Hooke’s law relates 1D normal stress to 1D extensional strain through two 
material parameters introduced previously: the modulus of elasticity E’, also called Young’s modulus 
and Poisson’s ratio v. The modulus of elasticity connects axial stress o to axial strain €: 


o-~Ee, E=-, €=—=. (5.3) 
€ E 


Poisson’s ratio v is defined as ratio of lateral strain to axial strain: 


lateral strain lateral strain 


(5.4) 


axial strain axial strain | 


The — sign is introduced for convenience so that v comes out positive. For structural materials 
v lies in the range 0.0 < v < 0.5. For most metals v © 0.25-0.35. For concrete and ceramics, 
v © 0.10. For cork v © 0. For rubber, v ~ 0.5 to 3 places. A material for which v = 0.5 is called 
incompressible. 
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$5.4 HOOKE’S LAW IN 1D 
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FIGURE 5.5. Speciment for determination of elastic modulus EF and Poisson’s ratio v in the linearly elastic 
response region of an isotropic material, using an uniaxial tension test. 


For distribution of shear stresses and 
strains over the cross section, cf. Lecture 7 


Vey 

0 

0 
at all points in the gaged region. Both the shear stress T,,= Ty as well 
as the shear strain Y= Yyx vary linearly as per distance from the 
cross section center (Lecture 7). They attain maximum values on the 
specimen surface. 
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FIGURE 5.6. Speciment for determination of shear modulus G in the linearly elastic response region of an 
isotropic material, using a torsion test. 


§5.4.2. Shear Modulus In 1D Stress State 


The shear modulus G connects a shear strain y to the corresponding shear stress T: 


T 
=Gy, G=-, = 55 
t y 7 Y=G (5.5) 


“Corresponding” means that if y = yxy, say, then tT = T,,, and similarly for the other shear 
components. The shear modulus is usually obtained from a torsion test. It turns out that the 
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1009019111 eee eee |||) 
«—— gaged length 


FIGURE 5.7. Speciment for determination of coefficient of thermal expansion a by heating a 
tension test specimen and allowing it to expand freely. 


3 material properties E, v and G for an elastic isotropic material are not independent, but are 
connected by the relations 
E 


G2 Posi sne 
21 +) eee v 


which means that if two of them are known by measurement, the third one can be obtained from 
the relations (5.6). In practice the three properties are often measured independently, and the 
(approximate) verification of (5.6) gives an idea of “how isotropic” the material is. 


§$5.4.3. Thermal Strains In 1D Stress State 


A temperature change of AT with respect to a base or reference level produces a thermal strain 
er =a AT, (5.7) 


in which @ is the coefficient of thermal dilatation, measured in 1/°F or 1/°C. This is typically 
positive: « > 0 and very small: a << 1, of order 10~° for most structural materials. To combine 
mechanical and thermal effects in 1D, the strains are superposed: 


e=enter= +a AT, (5.8) 


The last expression is valid if the material is linearly elastic and obeys the 1D Hooke’s law. 


Ox= 0 and €,,=0 at 
undeformed reference state, 
then bar is heated by AT 


E, a, Vv constant 
over bar 


FIGURE 5.8. Heated bar precluded from axial expansion. This bar will develop 
a compressive axial stress called a thermal stress. 


Example 5.1. The bar AB shown in Figure 5.8 is precluded from extending axially. It has elastic modulus 


E and coefficient of dilatation a > 0. The stress o is zero when the bar is at the reference temperature T,.;. 
Find which axial stress o develops if the temperature changes to T = T,-¢ + AT. 


Since the bar length cannot change, the combined axial strain must be zero: 


ex =€= E +a AT =0, (5.9) 
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§5.5 GENERALIZED HOOKE'S LAW IN 3D 


1 
vt 
a7 


(eal mii 9 
1 

NI 

AE SS hoe 


Final re 
after test (2) 


Final shape Oz Final shape 
after test (1) after test (3) 


FIGURE 5.9. Derivation of three-dimensional Generalized Hooke’s Law for normal stresses and strains. 
Three tension tests are assumed to be carried out along {x, y, z}, respectively, and strains superposed. 


Solving for o gives 
o =—EaAT’. (5.10) 


Since EF and a are positive, a rise in temperature, i.e., AT > 0, will produce a negative axial stress, and the bar 
will be in compression. This is an example of the so-called thermally induced stress or simply thermal stress. 
It is the reason behind the use of expansion joints in pavements, rails and bridges. The effect is important in 
orbiting vehicles such as satellites, which undergo extreme (and cyclical) temperature changes from full sun 
to Earth shade. 


§5.5. Generalized Hooke’s Law in 3D 
§5.5.1. Strain-To-Stress Relations 


We now generalize the foregoing equations to the three-dimensional case, still assuming that the 
material is elastic and isotropic. Condider a cube of material aligned with the axes {x, y, z}, as 
shown in Figure 5.9. Imagine that three“tension tests”, labeled (1), (2) and (3) respectively, are 
conducted along x, y and z, respectively. Pulling the material by applying o,, along x will produce 
normal strains 


(1) = Oxx (1) _ V Oxx (1) _ V Oxx 
Eo ae Evy = Ee”? €-, = ee (5.11) 
Next, pull the material by o,, along y to get the strains 
(2) _ Fyy Q) _  VOyy Q) _  VOyy 
Evy = EE? a= ee a €-7 => eae (5.12) 
Finally pull the material by o,, along z to get 
(3) _ Fz 3). 2 Vez Gy = Oz 
ey = =” SSG Ey =— EB (5.13) 
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In the general case the cube is subjected to combined normal stresses 0,,, Oyy and o,,. Since we 
assumed that the material is linearly elastic, the combined strains can be obtained by superposition 
of the foregoing results: 
_ 04 04,0) _ Sx _ YO%y VO _ | 
Ce GT eee eae E — E — E = Oxy — V Oyy — VOzz). 

SMC OD = OE 

E E E E 
VOxx VOyy Oz 1 
—— — — + = = — (-V 04x — VOyy + 02z). 

E E EE (—v ox. vy #2) 
The shear strains and stresses are connected by the shear modulus as 


() (2) (3) _ 
Eyy = Sy Tey 1 Sy = —VOxx + Oyy — Vv Gu) : (5.14) 


1 2 3 
a et + ee + ee os 


Tye Es Ties alles Tee << 1s 


i me Shama Ria a, TN ay er, — Vex = Yaz = = Ee (5.15) 
The three equations in (5.14), plus the three in (5.15), may be collectively expressed in matrix form 
as 
1 v v 
2 E -E -E 0 O O 4 
XX wv T _ D 0 0 0 XxX 
Eyy BE OE E Oyy 
Be Be 0 oO O 
Re as E EE ,; Pee (5.16) 
Yxy 0 0 0 G 0 O Txy 
Vy OQ: 80 HO <Q Bs SG) Ie 
Vex G 1 Tex 
0 0 0 0 O G 


§5.5.2. Stress-To-Strain Relations 


To get stresses if the strains are given, the most expedient method is to invert the matrix equation 
(5.16). This gives 


Cx E(L—v) | Ev Ep <0:“O--O5\fex, 
Oyy Ev E(l-—v) — Ev 0 0 O Eyy 
oy Ev Ev EQd-v) 0 0 O € 
oe 7 0 0 0 G 0 0 ms ' Onn 
Ue 0 0 0 0 G 0 Vyz 
Tg 0 0 0 0 0 G Vex 
Here E is an “effective” modulus modified by Poisson’s ratio: 
a (5.18) 
(1 — 2v)(1 + v) 
The six relations in (5.17) written out in long form are 
o Sopa IG she + V €yy + V Ezz | 
“(1 = 2v)(1 + v) os ? as 
E 
w= Toma) et OO er als (5.19) 
E 


= XxX ie ’ 
rae (1 —2v)(1 + v) est veg ts C=) ea| 


Try =Gyyy, Tye= Gy, Ty = GYzx. 
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The combination 


Cay = 5 (Oxx + Oyy + Oz) (5.20) 


is called the mean stress, or average stress. The negative of ogy is the pressure: p = —Oqy. 


The combination €, = €,, + €yy + €,; 1s called the volumetric strain, or dilatation. The negative 
of €, 1s known as the condensation. Both pressure and volumetric strain are invariants, that is, 
their value does not change if axes {x, y, z} are rotated. An important relation between pressure 
and volumetric strain can be obtained by adding the first three equations in (5.19), which upon 
simplification and accounting for (5.20) and p = —oz, relates pressure and volumetric strain as 


E 


Se eK ey. 21 
a =o) S 2") 


P = 
This coefficient K is called the bulk modulus. If Poisson’s ratio approaches +, which happens for 
near incompressible materials, K — oo. 
Remark 5.1. In the solid mechanics literature p is also defined (depending on author’s preferences) as 


D = Oq = + (Ors + Oyy + 0,,), which is the negative of the above one. If so, p = +K €,. The definition 
DP = —Oq, is the most common one in fluid mechanics. 


85.5.3. Thermal Effects in 3D 


To incorporate the effect of a temperature change AT with respect to a base or reference temperature, 
add a AT to the three normal strains in (5.14) 


1 

= a Oxy — VOyy — V Ozz) +a AT, 
1 

Sy = e (-v OF Oye Oz) +a AT, (5.22) 
1 

én = 5 (—v oxy — VOyy +0z,) +A AT. 


No change in the shear strain-stress relation is needed because if the material is linearly elastic and 
isotropic, a temperature change only produces normal strains. The stress-to-strain matrix relation 
(5.16) expands to 


1 v v 
we |e a ae 
xXx Vv Vv xx 
2 Te Se, pe OR OOM 1 
yy i oe 0 0 0 yy 
RE nde) ae. ae ae Bettie aa | (5.23) 
Vey 6 0 o £ © 6 |} 0 
Yyz CG 1 Tyz 0 
ae 0 0 0 0G ; oe 0 
0 0 0 0 0 @ 
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Inverting this relation provides the stress-strain relations that account for a temperature change: 


Orx E(1—v) s Ev Ev 0 0 O Exx 1 

Oyy Ev E(1— v) _ Ev 0 0 O yy 1 

Gigi > Ev Ev E(l-v) 0 0 O é., | EaAT] 1 

Paes ¥ 0 0 0 G 0 0 Yry f= 2v: |0:)” 

Tyz 0 0 0 0 G 0 Vyz 0 

i 0 0 0 0 0 G Vx 0 
(5.24) 


in which E is defined in (5.18). Note that if all mechanical normal strains €,,, €,), and €,, vanish, 
but AT + 0, the normal stresses given by (5.24) are nonzero. Those are called initial thermal 
stresses, and are important in engineering systems exposed to large temperature variations, such as 
rails, turbine engines, satellites or reentry vehicles. 


§5.6. Generalized Hooke’s Law in 2D 


Two specializations of the foregoing 3D equations to two dimensions are of interest in the appli- 
cations: plane stress and plane strain. Plane stress is more important in Aerospace structures, 
which tend to be thin, so in this course more attention is given to that case. Both specializations 
are reviewed next. 


§5.6.1. Plane Stress 


In this case all stress components with a z component are assumed to vanish. For a linearly elastic 
isotropic material, the strain and stress matrices take on the form 


Exx Yay 0 Oxx Txy 0 
VYyx yy 0 | ’ Tyx  Oyy 7 (5.25) 
Oi 0. Sey 0. 0° 0 


Note that the €,, strain, often called the transverse strain or thickness strain in applications, in 
general will be nonzero because of Poisson’s ratio effect. The strain-stress equations are easily 
obtained by going to (5.14) and (5.15) and setting o,, = Ty, = T;x = 0. This gives 


; v 
tases E (CCS Voy). eg = E (Vox toy), ee = —E (xx + Oy) (5.26) 
Yxy = o Ne ta = 0: 


The matrix form, omitting known zero components, is 


1 v 
ey Be aS ihe 
Vv XX 
€ —7 7 O 
ce =|_ a a i Gy | ; (5.27) 
zz E E Ly 
ee 0 0 
Inverting the matrix composed by the first, second and fourth rows of the above relation gives the 
stress-strain equations 23 - 
Cs E Ev 0 Exx 
a = [és E | fs : (5.28) 
Txy 0 0 G Yay 
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in which E = E/(1 — v?). Written in long hand, 


E 
To & Meyer). Tey = G Vee (5.29) 


Ox = To Se ale Eyy), Oyy = 


1 


If temperature changes AT are considered, the foregoing equations acquire extra terms on the right: 


1 v 
ex Ee 2 I 
: 22¥~ IS — 4G Ox 1 
yy | E FE. Oyy | +a AT ; (5.30) 
Ezz —7 —p- O | 
Very | 2 = 1 | ay 0 
0 0 G 
and ~ ~ 
Ors BEY OT Ge) oar 2 
Oy |=] Ev E O]] ey ]- Lil (6.31) 
l-—v 
Ly 0 0 Gilly 0 


$5.6.2. Plane Strain 


In this case all strain components with a z component are assumed to vanish. For a linearly elastic 
isotropic material, the strain and stress matrices take on the form 


Exx YVry 0 Oxx Try 0 
Vox €yy 7 : ye Oy 0 | (5.32) 
0 0 0 0. yes 


Note that the o,, stress, which is called the transverse stress in applications, in general will not 
vanish. The strain-to-stress relations can be easily obtained by setting €,, = yz = yz = Oin 
(5.19). This gives 


E 

Oxx = (ony a) [d — V) Exe + Vey | ; 
E 

De (Oey) ential: (5.33) 
E 


Ozz = al = 2v)(1 a v) [v Exx + V yy | ’ 


tLe OVs6 Gy =U. = 0. 


which in matrix form, with the zero components removed, is 


Or E dL —v) . Ev 0 ‘ 

ora Ev E(l—v) 0 ies 

Ge || - Ev Ev 0 ay: |? Oey 
Try 0 0 Gir 


Inverting the system provided by extracting the first, second and fourth rows of (5.34) gives the 
stress-to-strain relations, which are omitted for simplicity. 


The effect of temperature changes can be incorporated without any difficulty. 
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Backup Material Only — Not Covered In Lecture 


Initial (reference) 
diameter Dp, under 
inflation pressure pj 


oe 


For simplicity assume 
aspherical ballon shape 
for all calculations 


inflation 
pressure 


P; = Pot P 


Final (deformed) 
diameter D, = B D, 


| ee B Dy) ———> 


FIGURE 5.10. Inflating balloon example problem. 


§5.7.. Example: An Inflating Balloon 


This is a generalization of Problem 3 of Recitation #2. The main change is that all data is expressed 
and kept in variable form until the problem is solved. Specific numbers are plugged in at the end. 


This will be the only problem in the course where some features of nonlinear mechanics appear. 
These come in by writing the governing equations in both the initial and final geometries, without 
linearization. 


$5.7.1. Strains and Stresses in Balloon Wall 


The problem is depicted in Figure 5.10. A spherical rubber balloon has initial diameter Dp under 
inflation pressure po. This is called the initial configuration. The pressure is increased by p so 
that the final pressure is py = po + p. The balloon assumes a spherical shape with final diameter 
Dy = B Do, inwhich B > 1. This will be called the final configuration. The initial wall thickness is 
to << Dp and the final thickness is ty. Since the balloon geometry is assumed to remain spherical 
for simplicity, we can apply to both configurations the stress formulas for the thin-wall spherical 
vessel derived in Lecture 3. 


The strains (but not stresses) are assumed to be zero in the initial configuration. The average 
circumferential extensional strain assumed in the final configuration depends on whether we take 
the Lagrangian or the Eulerian strain measure, which are designated by ¢“, and €Z,, respectively. 


Obviously 


(Df — Do) = 


L 
Er = 
f ma Do 


pl, 6 =tee ey = (5.35) 


§5.7 EXAMPLE: AN INFLATING BALLOON 


Since the balloon is assumed to remain spherical and its thickness is very small compared to its 
diameter, the above strains hold at all points of the balloon wall, and are the same in any direction 
tangent to the sphere. If we choose the sphere normal as local z axis, the wall is in a plane stress 
state. 


Next we introduce material laws. We will assume that rubber obeys the two-dimensional, plane 

stress generalized Hooke’s law (5.31) with respect to the Eulerian strain measure, with effective 

modulus of elasticity E and Poisson’s ratio v! Setting €xx = €y = ef and y,, = 0 therein and 

accounting for the initial stress 09, we obtain the inplane normal stress in the final configuration: 
fp eh Ge) ope EE, BG AG) 

OS Vyy = Or = oo Sse ver) = 09 + —— €; = 09 + —— ——. : 
- ed i aie ee f Lay v 
The normal inplane wall stress is the same in all directions, so it is called simply oo and o,, for 
initial and final configurations, respectively. The inplane shear stress t vanishes in all directions. 


Assume Dp, fo, E and v are given as data. An interesting question: what is the relation between p 
(the excess or gage pressure) and the diameter D¢ = BDo? And, is there a maximum pressure that 
will cause the balloon to burst? 


§5.7.2. When Will the Balloon Burst? 


To relate p and # it is necessary to express the wall stresses oo and of in terms of geometry and 
internal pressure. This is provided by equation (3.10) in Lecture 3, derived for a thin-wall spherical 
vessel. In that equation replace p, R and t by quantities in the initial and final configurations: 


on, — PoRo _ PoPo e _ prRe — (Potp)Dye _ (P0+ P)B Do 
a 2to Ato : f 2 At Ate . 


(5.37) 


All quantities in the above expressions are known in terms of the data, except ty. A kinematic 
analysis beyond the scope of this course shows that 


th = c +2v (= = i) to. (5.38) 


We can check (5.38) by inserting two limit values of Poisson’s ratio: 

v=0: tr = ty. This is correct since the thickness does not change. 

v=1/2: tr = t)/B?. Is this correct? If v = 1/2 the material is incompressible and does not 
change volume. The initial and final volume of the thin-wall spherical balloon are 
Vo = 1 Dp to and Vy = 1D ty = 1B? Doty, respectively. On setting Vo = V> and 
solving for t we get tp = to/B. 

To obtain p in terms of £, replace (5.38) into (5.37), equate this to (5.36) and solve for p. The 

result provided by Mathematica is 


a= 4Eto(1 — B)(2v + B*(1 — 2v)) + DopoBU — v)(4v + 672 — B — 4v)) 


DoBid — v) (5.39) 


' This is a very rough approximation since constitutive equations for rubber (and polymers in general) are highly nonlinear. 
But getting closer to reality would take us into the realm of nonlinear elasticity, which is a graduate-level topic. 
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Lecture 5: STRESS-STRAIN MATERIAL LAWS 


Fi v=0 


4 B=D;/Do 


is 2 2 3 35 


FIGURE 5.11. Inflating pressure (in MPa) versus diameter expansion ratio 
B = Dy /Do foraballoon with E = 1900 MPa, Dp = 50mm, po = 0 MPa, 
1 


to = 0.18 mm, 1 < 6 < 4 and Poisson’s ratios v = 0 and v = ie 


This expression simplifies considerably in the two Poisson’s ratio limits: 


AE ty (B — 1) + Do po (2 — 
eee to (B a = Po(2— B)B (5.40) 


8E to (8B — 1) + Do po (2 — B°) B 
Ply=ip2 = Do B4 


(5.41) 


Pressure versus diameter ratio curves given by (5.40) and (5.41) are plotted in Figure 5.11 for 
the numerical values indicated there. Those values correspond to the data used in Problem 3 of 
Recitation 2, in which v = 1/2 was specified from the start. 


Rubber (and, in general, polymer materials) are nearly incompressible; for example v ~ 0.4995 for 
rubber. Consequently, the response depicted in Figure 5.11 for v = 1/2 is more physically relevant 
than the other one. 

Do the response plots in Figure 5.11 tell you when an inflating balloon is about to collapse? Yes. 
This is the matter of a (optional) Homework Exercise. 
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